Structure and dynamics in glass-formers: predictability at large length scales 
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Dynamic heterogeneity in glass- formers has been related to their static structure using the con- 
cept of dynamic propensity. We re-examine this relationship by analyzing dynamical fluctuations 
in two atomistic glass-formers and two theoretical models. We introduce quantitative statistical 
indicators which show that the dynamics of individual particles cannot be predicted on the basis 
of the propensity, nor by any structural indicator. However, the spatial structure of the propensity 
field does have predictive power for the spatial correlations associated with dynamic heterogeneity. 
Our results suggest that the quest for a connection between static and dynamic properties of glass- 
formers at the particle level is vain, but they demonstrate that such connection does exist on larger 
length scales. 

PACS numbers: 05.10.-a, 05.20.Jj, 64.70.Pf 



The future ain't what it used to be — Y. Berra 



I. INTRODUCTION 

Supercooled liquids near the glass transition have 
paradoxical physical properties [1, 2]. They exhibit a 
range of peculiar dynamical features that have been 
attributed to spatially het erog eneous dynamical relax- 
ation 0, H, i, U i, 0, i, i, Ea but their structure, 
as measured by two-point correlation functions, appears 
homogeneous and unspectacular. Theoretical pictures of 
the glass transition assume different kinds of connections 
between static and dynamic properties. For example, in a 
picture based on dynamical facilitation [H, [lH , one postu- 
lates the existence of mobile and immobile regions, with 
the implicit assumption that these regions have a struc- 
tural origin. Frustration-based theories [13] infer dynam- 
ical behavior by assuming the existence of domains with 
a preferred local order. Alternatively, one can attempt 
to connect static and dynamical properties through the 
configurational entropy [3, [lH , through two-point den- 
sity correlations [I^, through elastic properties [I?, 
or through the idea of a rough energy landscape [7, 
The extent to which these connections can be objectively 
established in experiments and computer simulations is 
an important criterion for evaluating different theoretical 
pictures. 

To address this point, Harrowell and co-workers intro- 
duced the isoconfigurational ensemble [20] , which isolates 
the effect of the liquid structure on its dynamical fluctu- 
ations. In this statistical ensemble, dynamical observ- 
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ables known as 'propensities' are obtained by averaging 
over independent trajectories from the same initial con- 
figuration. Fixing the initial particle positions preserves 
structural information on all length scales, and allows the 
influence of structure to be separated from intrinsically 
dynamical fluctuations. The difficult task of connecting 
structure to dynamics is then broken into two seemingly 
simpler ones: first connect structure to propensity, then 
propensity to dynamics. 

This idea has been exploited in several works [lO, [2l|, 
S S S m, S [13, which focussed primarily on the 
first stage of the problem, and aimed at finding corre- 
lations between propensity and local or non-local struc- 
tural quantities. In this paper, we concentrate instead 
on the second stage, the connection between propensity 
and dynamics, which has received relatively little atten- 
tion. After all, if propensity and dynamics only had 
weak connections, the isoconfigurational ensemble would 
not be such a useful tool. The isoconfigurational aver- 
age preserves all fluctuations that have any connection 
to the initial structure, but fluctuations whose origin is 
inherently dynamical are averaged away. Thus, propen- 
sities capture the physically relevant dynamical fluctua- 
tions if and only if these fluctuations are structural in 
origin [20], in which case dynamics can be predicted 
from knowledge of the structure. In concentrating on 
the connection between structure and propensity, previ- 
ous works m, m, m, m, m, [23 have assumed that this 
is indeed the case. 

Here, we analyze quantitatively the implicit assump- 
tion of predictability. We establish that the dominant 
single-particle dynamical fluctuations are intrinsically 
dynamical, and not linked with liquid structure. How- 
ever, we find that the sizes and shapes of mobile and im- 
mobile regions can be predicted from the propensity, and 
these collective dynamical fluctuations do indeed have a 
structural origin. Although the relevance of intrinsically 
dynamical fluctuations has been discussed in a qualita- 
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tive way by Widmer- Cooper and Harrowell [10, [2l|, the 
subtle length scale dependence of the correlation between 
structure and dynamics, and therefore of predictability, 
was not discussed nor anticipated in previous work. 

After defining our models in Sec. [TTl we discuss single 
particle and collective dynamics in Sees. IIIII and [IV| com- 
bining illustrative snapshots with quantitative analysis, 
and comparing our results with the behavior in schematic 
models. We conclude in Sec. [Vl and identify directions 
for future study. 



II. MODELS 

We present numerical data for two atomistic glass- 
formers: a Lennard- Jones (LJ) binary mixture [28], and 
a model of silica due to Beest, Kramer and van San- 
ten (BKS) [1^. We use Monte Carlo dynamics, which 
have been shown to yield a dynamical behavior in excel- 
lent agreement with Newtonian dynamics when time is 
scaled appropriately [sO, [3l| . We measure time in Monte 
Carlo sweeps, and other units are as in Refs. [sO, Isij . 
The LJ system has 1000 particles, and we study temper- 
atures between 1.0, where the system is a simple liquid, 
and 0.47, where the relaxation time has increased by a 
factor of approximately 1000, and the system is in the 
glassy regime. For reference the mode-coupling temper- 
ature for this system is = 0.435. The BKS system has 
1008 atoms, and the relaxation time also spans around 
three decades, from a liquid state around 6000 K to a 
glassy one around 3000 K. 

We also study the one-spin facilitated Fredrickson- 
Andersen (1-FA) model p^], using Monte Carlo simu- 
lations. It represents a simple model of a dynamically 
heterogeneous material, with a few mobile regions that 
facilitate motion in immobile regions nearby [5|] . Defining 
spins Hi G {0, 1} on the site of a lattice, we identify sites 
with = 1 as 'mobile' and those with = as 'immo- 
bile'. The dynamics of the system obey detailed balance 
with respect to a trivial energy function E = Ui. The 
non-trivial behavior of the model comes from a dynam- 
ical constraint: spins are allowed to fiip only if at least 
one of their neighboring sites is mobile. At low tempera- 
tures, T < 1, and low spatial dimension, d < 2, the model 
displays glassy features such as transport decoupling [32[ 
and an increasing dynamical length scale 

We use Ci{t) to denote a general dynamical object at- 
tached to particle such as fi{t) = cos{k ■ [ri{t) — Ti(0)]) 
or ^i{t) = \n{t) - n(0)|. For fi{t), we use \k\ = 6.7 for 
the LJ model and \k\ = 1.8 A~^ for silica. These wave 
vectors correspond to the location of the first diffraction 
peak in the LJ system, and to the pre-peak of the struc- 
ture factor of silica. The quantity fi{t) is therefore a 
measure of local relaxation in these liquids. In lattice 
models, we consider Pi(t), the persistence function on 
site z, which takes the value of unity if spin i has not 
fiipped in the interval [0,t], and zero otherwise. 

We define an isoconfigurational average (• • •)iso in 



which the initial positions of all particles are held 
fixed [13], but dynamical trajectories are made indepen- 
dent through the use of different random numbers in the 
Monte Carlo trajectories [25[. The dynamic propensity 
is then defined by (Ci(t))iso HO]- Equilibrium ensemble 
averages are denoted by E[- • •]. As usual, we define the 
structural relaxation time by E[/^(rQ,)] = 1/e. 



III. SINGLE PARTICLE DYNAMICS 
A. Atomistic systems 

We begin with a qualitative analysis of fiuctuations 
within the isoconfigurational ensemble. Fig. [1] shows a 
scatter plot for the LJ system, comparing squared parti- 
cle displacements in a single run, /if (t), with their corre- 
sponding propensities, (/i|(t))iso. Points near the dashed 
line represent particles whose displacements in this single 
run are close to their isoconfigurational averages. Given 
the large scatter seen in Fig.[TJ it is clear that the propen- 
sity cannot be used to predict the actual value of the 
displacement in a single run. That is, even when all as- 
pects of the initial structure are held constant, there are 
large fiuctuations in the single-particle dynamics, which 
appear to be more important than particle-to-particle 
fiuctuations of the propensity. These large dynamical 
fiuctuations mean that the propensities have very little 
predictive power for the single-particle dynamics. 

Fig. [U also shows the distribution of particle (loga- 
rithmic) displacements, P(log/i), and of propensities, 
P(log(/i)iso). The former, which is directly related to the 




FIG. 1: Data for the LJ system in the glassy regime (T = 
0.47). Inset: Scatter plot, comparing squared displacements 
in a single run, /if(Ta), with their average over runs from 
the same initial condition, (//^(tq, ))iso- Note that the lack 
of correlation between propensity and dynamics necessitates 
the use of different scales on the two axes. The dotted line 
is /if (tq;) = {lJii{T a)) ISO- Main: The broad distribution of par- 
ticle (logarithmic) displacements, P(log /x) (full line) arises 
from the coexistence of mobile and immobile particles, illus- 
trated with dashed lines (Gaussian distributions). This is 
compared with the much narrower distribution of propensi- 
ties, P(log(/i)iso). 
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van-Hove function, has a broad, non-Gaussian shape, re- 
flecting the coexistence of mobile and immobile particles 
in the liquid 0, [ll|, as illustrated by fltting the large- /i 
and small- /i parts of the distribution with two distinct 
Gaussian distributions. The latter distribution (particle 
propensities) is much narrower and structureless. In par- 
ticular, the distinction between fast and slow particles 
is no longer apparent. This most distinctive feature of 
dynamic heterogeneity [Ij] is therefore not structural in 
nature, but is intrinsically dynamical. 

The details of the distribution of the propensity shown 
in Fig. [T] are different from those shown in [20|,|2l|. We at- 
tribute this to differences between our model systems (in 
particular their different dimensionalities and tendencies 
to crystallize) , and to our use of slightly different observ- 
ables (we use simple distances while squared distances 
were used in [lO, [2l|, emphasizing particles with large 
displacement). For our purposes, the important feature 
is that the distribution of propensities is much narrower 
than that of bare displacements, which is consistent with 
earlier results [33^ . 

We now support these qualitative statements by quan- 
titative measures. In order to disentangle structural and 
dynamical sources of fluctuations, we deflne three vari- 
ances: 



Aj?°(t) 
Ac(t) 



E[{c!m 



' 1 

isoj 



-E^ [CM, 



(1) 



SO that Ac{t) = Ag°(t) + Sc{t); we use the short-hand 
notation E^[- • •] = (E[- • -J)^ and the subscript indicates 
the dynamic observable of interest. Thus, Sc{t) measures 
particle to particle fluctuations of the propensity [l^ and 
captures the structural component of the fluctuations. 
Ag°(t) measures the fluctuations of Ci{t) between differ- 
ent runs at flxed initial conflguration and captures there- 
fore the dynamical component of the fluctuations. Their 
sum Ac{t) naturally measures the total amount of fluc- 
tuations. 

These quantities are deflned in the same spirit as the 
ensemble-dependent susceptibilities of Refs. [10|. To as- 
sess the influence of a given variable (the structure) 
on the fluctuations measured by Ac(t), we use a con- 
strained statistical ensemble (the isoconflgurational en- 
semble) where this variable is kept flxed, and measure 
a 'restricted' variance, Ag°(t). The difference between 
restricted and unrestricted variances is fc(^), which ac- 
counts for the fluctuations of the restricted variable [TO] • 
In our case, the relative sizes of 6c {t) and Ag°(t) quan- 
tify the relative influence of structural and dynamical 
fluctuations. It is natural to introduce the dimensionless 
ratio 



Rc{t) 



ML 



(2) 



with < Rc{t) < 1. Small values of Rc{t) mean that the 
structure has little effect on the single particle dynamics. 



0.4 
0.3 
0.2 
0.1 



1 ' ' ' 




1 




LJ - -0 - 
BKS - Si 
BKS - ---A-- 














III. 1 


.-..-.4-O-0 ^/(^) 
1 



102 



10^ 



10^ 

Try. 



10^ 



FIG. 2: Ratios Rf{Ta) for the LJ and silica systems, and 
R4,f{ra) for the LJ system. Temperatures are 0.47 < T < 
1.0 (LJ) and 3000K < T < 6100K (BKS). Rf{t) is small 
and grows slowly with decreasing T, while i?4,/(t) is large 
and grows steadily. This suggests that connections between 
structure and dynamics are only significant on large length 
scales. 



large values mean instead that structure is a very good 
predictor of the dynamics. 

We present results for the ratio Rf{ra) for LJ and BKS 
glass- formers in Fig.O This ratio vanishes at high T, and 
grows slowly as increases, reaching about 4% in both 
systems at the lowest temperature studied (when has 
grown by about 3 decades). The smallness of Rf{ra) 
conflrms the impression gained from Fig. [TJ structure at 
time has little influence on the dynamics of individual 
particles at times r^. 



B. Comparison with a schematic model 

It is instructive to evaluate Rcif) in a simple model 
of a dynamically heterogeneous liquid. We suppose that 
particles diffuse independently with diffusion constants 
Di that depend on their initial environments. We also 
assume that effects of the initial environment decay on 
timescales comparable with Ta- Then, each particle has 
m))iso = e-^-'^'*, and 



F{t)=E[m] 



dAi/(A)exp(-At), (3) 



where = Dik^ is a rescaled diffusion constant dis- 
tributed according to ^(A). Using the definition of /i(t), 
we have fi{tf = (l/2){cos(2fe • [ri{t) - ri{0)]) - 1}, and 
we find 



5f{t) = F{2t)-F\t), 
Afit) = (l/2)[l + F(4t)-2F2(t)]. 



(4) 



In a dynamically homogeneous system, all particles 
have the same diffusion constant, F{t) decays exponen- 
tially, and Sf{t) = = Rfit). On the other hand, 
consider a heterogeneous system with equal populations 
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of fast and slow particles, whose rescaled diffusion con- 
stants are Ai and A2. Using Eq. ([3|), this leads to a two- 
step decay: F{t) = (1/2) [g-^^^ + e'^^*]. In this case, 
Rf{t) has a non-monotonic time dependence, vanishing 
at small and long times, with a maximal value near 50% 
during the plateau of F{t) (that is, for times such that 
Xi^ <C t <C A^^). This indicates that correlations be- 
tween structure and dynamics are generically maximal 
during plateaux of F{t) [23]; our atomistic simulations 
are also consistent with this indication. 

Our analysis of this highly schematic model demon- 
strates that large values of Rf{t) are obtained in the 
presence of reproducibly fast and slow particles. Compar- 
ing this reference theory with the small values for Rf{t) 
shown in Fig. [21 it follows that the correlation between 
structure and single-particle dynamics is weak in atom- 
istic systems and does not seem to dramatically increase 
when temperature gets smaller. 

We conclude that the search for a connection be- 
tween static and dynamic properties at the single-particle 
level is in vain. Certainl y, n o such connection has been 

found [23, [m, [13, [13, His, m, [13 . 

IV. COLLECTIVE DYNAMICS 

Having ruled out the structural origin of one impor- 
tant aspect of dynamic heterogeneity, we now address a 
different question: Do spatial fluctuations of the propen- 
sity carry meaningful information on the geometry and 
spatial extent of the dynamic heterogeneities in super- 
cooled liquids? We will show that while it is not possible 
to use the structure to predict whether a given particle 
will be fast or slow in a single run, it is possible to tell if 
it belongs to a fast or slow region. 



A. LJ system 

We begin by generalizing the variances in ([1]) to those 
of global dynamic quantities, Cg{t) = N~'^ Ci{t). We 
define: 

54,c(t) = iV{E[(Cg(t))L] -E2[Cg(i)]}, 
xZcit) = A^{E[{C|W)i.o-(Cg(t))L]}, 
X4.c(t) = N{E [(C|(t))i,o] - [C,m, (5) 

with X4,c(^) = ^4,c(0 + X4 cl^)- The usual four-point 
susceptibility X4,c(^) measures the size of collective dy- 
namical fluctuations [1, By analogy with Sc{t) and 
AJ3°(t), we identify S^^cit) as the contribution to X4,c(^) 
associated with the structure, and X4 c(^) intrin- 
sically dynamical contribution. Then, in analogy with 

Rc{t), 

R^Mt) = (6) 



is a dimensionless measure of the effect of the structure 
on collective aspects of the dynamics. In Fig. [21 we show 
that the contribution of ^4,/(ra) to the dynamic fluctua- 
tions is about 35% at T = 0.47 (to be compared to the 
4% found for single particle heterogeneity). Moreover, 
RAjijoc) grows steadily when T decreases. This quan- 
titative measurement shows that when the system is in 
the glassy regime, collective dynamical fluctuations are 
indeed quite reproducible in repeated runs from the same 
initial configuration. We conclude that connections be- 
tween structural and dynamical properties are significant 
on these length scales, and that the connection even gets 
stronger as the relaxation time increases. 

To investigate the fluctuations associated with this ef- 
fect, we interpolate between the single-particle and col- 
lective dynamics, by averaging over a length scale t 

E;l(|.-' -r,|) ' (7) 

By coarse-graining in this way, we can investigate how 
the reproducibility of dynamical fluctuations varies with 
length scale. 

In Fig. [31 we use snapshots of the system to illustrate 
that large scale features of the dynamics are indeed re- 
producible from run to run. For a representative ini- 
tial configuration at a low temperature. Fig. [3^ shows 
the LJ particles with the largest values of the propensity 
{fi{^a)) ISO (these are the particles which are slow on av- 
erage). In all panels of Fig. [3l we show about 1/3 of the 
particles: this threshold is small enough to give clear im- 
ages, but large enough to avoid placing undue emphasis 
on rare fluctuations. 

Figs. [3)3-d show particles with the largest values of 
fi{t,£ = 1) in three individual runs. Loosely speaking, 
the coarse-graining scale i = 1 means that measures 
how much motion is associated with a particle and its 
nearest neighbors, on the time scale t (the choice of this 
length scale is discussed below). Thus, Figs. [Sb-d show 
particles that are located in a relatively immobile envi- 
ronment during these three trajectories. Similar clusters 
of immobile particles are apparent in all three trajecto- 
ries, and these clusters correlate quite well with the clus- 
ter of slow particles that is observed in the propensity 
map in Fig. [3^1. Since the slow behavior of these clusters 
is reproducible in independent runs from the same initial 
condition, it surely must have a structural origin. This 
is consistent with the rather strong coupling of structure 
and collective dynamics that was identified in Fig. [2l We 
also note in passing that the coarse-grained pictures in 
Figs. [3)3-d are computationally much cheaper than cal- 
culating the propensity field. 

As discussed above, the coupling between structure 
and dynamics at the single-particle level is weak. This is 
further illustrated in Fig. [3^, where we have not coarse- 
grained, but simply identified slow particles by their val- 
ues of /i(rc^), again using the same initial condition. In 
this case the immobile cluster that is apparent in Fig.[3^- 
d is obscured by large intrinsically dynamical fluctua- 
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FIG. 3: a: Image of the 300 particles with the smallest propen- 
sity for motion (largest values of {fi{ra))iso), for a represen- 
tative initial condition in the LJ system in the glassy regime 
(T = 0.47). b-d: Images of the 300 particles with largest val- 
ues of the coarse-grained observable fi{Ta,i = 1), in three 
representative runs from the same initial condition. Spatial 
structure similar to that of (a) is found in all cases, e: As (c), 
without coarse- graining (i = 0). f: Image of the 300 parti- 
cles with the smallest values of the coarse-grained potential 
energy, ei{i = 2), for the same configuration. The resulting 
structure is correlated with that of (a). 



tions. Comparing Fig. [3^ with Figs. [3)3-d shows that the 
effect of coarse-graining on the short length scale i = 1 
is quite effective in suppressing these fluctuations, allow- 
ing the slow cluster to become apparent. When coarse- 
graining in this way, we must also ensure that the length 
scale i is smaller than the correlation length associated 
with the dynamically correlated clusters, or else these 
clusters will themselves be obscured. We establish below 
(Fig. [4]) that the dynamically correlated clusters have a 
length scale ^p^°p :^ 2 at this temperature. Thus, while 
it would be desirable to have well-separated length scales 
£ and ^p^°p, we can at least establish that i < ^p^°p, 
as required for the consistency of our analysis. As ex- 
pected, we find that on further increasing the coarse- 
graining scale the structure of the immobile clusters 
in Figs. [3)3-d is still apparent, but the ability to resolve 




FIG. 4: Top: Images of the 300 particles with the largest val- 
ues of (/z(rc))iso at T = 1.0, 0.6, and 0.47 (left to right) in the 
LJ system. Increasing clustering of propensity fluctuations 
is evident. Bottom: Structure factor of propensity fluctua- 
tions at several T. The inset shows the extracted lengthscale 
^^^°^(T), which grows as the temperature decreases. 



their shape is lost. 

Returning to the spatial correlations of the propen- 
sity, we interpret S4j{t) as a dynamic susceptibility as- 
sociated with spatial fluctuations of the propensity, b; 
analogy with the four-point susceptibility X4,/(^) 
That is, defining the fluctuations of the propensity by 
{^fi{t))iso = {fi{t))iso — F{t)^ then the spatial correlation 
function of the propensity is 

x,5(n(0)-r,(0)-r)]. (8) 

The associated structure factor is the Fourier transform 
of this function: 



^prop 



dr e 



(9) 



and the associated dynamical susceptibility is ^4,c(^) = 
5'4^°^(g ^ 0, t). By analogy with the four-point suscepti- 
bility, we expect (54,c(^) to be proportional to the number 
of particles associated with collective fluctuations of the 
propensity. 

Like XA^(f)i the susceptibility of the propensity S/^^ci^) 
is a non-monotonic function of time that peaks near Tq,. 
We find that the height of this peak grows as T decreases, 
suggesting increasing spatial correlations of the propen- 
sity. This is illustrated by the images in Fig. [H which 
show that the particles with smallest propensities are in- 
creasingly spatially clustered as T decreases. To confirm 
this visual impression, we present in Fig. [Hour numerical 
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data for S^^'^^{q^ t) measured at different temperatures in 
the LJ system. The structure factor of the propensity is 
similar to the structure factor of dynamic heterogeneity, 
with the striking appearance of a smah-g peak. Fohow- 
ing earlier work we estimated [34] the correlation length 
of propensity fluctuations ^p^op(T), as shown in the inset 
of Fig. m It has a clear, but rather slow, growth when 
T decreases, compatible with determinations of a typical 
lengthscale of dynamic heterogeneity [H, 0, [T^ . 

These measurements confirm that spatial correlations 
of the propensity increase when temperature decreases, 
just as the length scale of dynamic heterogeneity does. 
Therefore, we find the intriguing result that the value of 
the propensity of any individual particle is only weakly 
correlated to its dynamical behavior, but the spatial cor- 
relations of these propensities do carry information about 
spatially heterogeneous dynamics. In short, the spatial 
structure of the propensity maps in Ref. ^ is important, 
but the color assigned to any specific particle is not. 



B. Comparison with schematic model 

We now generalize the simple model of Sec. lIIIBl to in- 
clude spatial correlations. Following [4], we assume that 
particles diffuse independently, but with diffusion con- 
stants that are correlated over large spatial regions, each 
containing nc ^ 1 particles. Since particles diffuse in- 
dependently in any given run of the dynamics, it follows 
that {fi{t)fj{t))iso = {fi{t))iso{fj{t))iso for j 7^ i, and 
hence: 

S^j{t) = n,[F{2t)-F\t)] 

Xr,}(t) = (l/2)[l + F(4t)-2F(2t)]. (10) 

In this simple model, both d/^jit) and X4,/(^) scale with 
nc, while the isoconfigurational susceptibility X4 /(^) is 
not sensitive to spatial correlations of the mobility and 
remains (9(1). Thus, X4,/(^) ^ ^^ji^) > xij(^), and 
hence R/^j :^ 1. This model shows that R4j{t) becomes 
large if the lengthscale for dynamic heterogeneity is pri- 
marily structural in nature. The results of Fig. [2] there- 
fore indicate that this is the case for the LJ system. 



C. Comparison with a kinetically constrained 
model 

We end with a brief discussion of the 1-FA model, 
which gives a useful theoretical insight into the quanti- 
ties discussed above. As recalled in Sec. [Ill the model de- 
scribes a dynamically heterogeneous material in which a 
few mobile excitations diffuse through an immobile back- 
ground. Our simulations indicate that the ratios Rp{t) 
and R4^p{t) have limiting forms at low T, which depend 
on dimensionality, d. In Id, the structure is very strongly 
correlated with the dynamics both locally and globally: 
[Rp{r)^ R4^p{r)] ^ [0.5,0.7]. However, in 3d^ the single 



site ratio vanishes, Rp{r) ^ 0, while the global ratio is 
quite large, i?4,p(r) ~ 0.4. This occurs because the set 
of sites visited by a given excitation in 3d varies enor- 
mously from run-to-run. Using the initial positions of 
excitations to predict which sites will relax first in a given 
run is impossible. However, collective observables reveal 
that the rate of relaxation is reproducibly enhanced in 
regions with relatively many excitations (see also [25l]). 

This decoupling between the local and global ratios 
illustrates a situation in which the relation between dy- 
namics and structure is only statistically significant at 
large length scales. Interestingly, this result is somewhat 
similar to that shown in Fig. [2] for the LJ system, al- 
though the microscopic mechanisms at work are presum- 
ably different. Moreover, preliminary studies indicate 
that the strong length scale dependence of predictabil- 
ity found in the 1-FA model is much less pronounced in 
models where kinetic constraints are stronger. In these 
other models, it appears that the dynamics on all length 
scales is strongly constrained by the initial structure. 



V. OUTLOOK 

We have investigated the degree to which structural 
fluctuations in glass-forming liquids influence their dy- 
namical fluctuations. We defined R^j and Rf which are 
quantitative measures of this influence, on long and short 
length scales respectively: they differ by nearly an order 
of magnitude in the LJ system at the lowest temperature 
considered. Thus, the influence of structure on dynamics 
is much stronger on long length scales than on short ones. 
This unexpected finding constitutes our main result. 

Our work does not reveal which structural features are 
responsible for dynamic heterogeneity, but they do show 
that the search for such an observable should be under- 
taken at a coarse-grained level. This is consistent with 
recent studies [13, lli, H, [25[ . 

As a first step towards identifying a suitable coarse- 
grained structural quantity, we exploit the fact 0, 
that potential energy is correlated with dynamical het- 
erogeneities (although correlations of the energy remain 
short-ranged). We compare fluctuations of the coarse- 
grained energy field, ei(^), with those of the propen- 
sity. By coarse-graining on a length scale £ = 2, so that 
£ ^ ^prop^ average away local fluctuations. We ob- 
tain a field that varies in space on a similar length scale 
to that of the dynamical propensity, and which reflects 
the average energy of different regions of the system. In- 
terestingly, we find that regions with small energy are 
correlated with regions of low propensity for motion, as 
shown in Fig. [3^. While a more quantitative analysis 
is necessary before drawing firm conclusions, this corre- 
lation between energy and propensity is consistent with 
the strong local correlations between energy and dynam- 
ics demonstrated in 0, [13] • As an alternative to the 
energy, another promising route to a connection between 
structure and dynamics is provided by the presence of 
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extended modes characterizing the vibrational spectrum 
of amorphous materials [13, [Si? it would also be 
interesting to study connection between propensity and 
the locally ordered regions discussed in Ref. [35]. 

In any case, identifying the non-local structural fea- 
tures that are associated with mobile or immobile regions 
of glass- formers remains a central challenge. 
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